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Abstract 

Bound-state-like wave functions are used to determine the scattering matrix corresponding to low 
energy N — d and p— 3 He collisions. To this end, the coupled channel form of the integral relations 
derived from the Kohn variational principle is used. The construction of degenerate bound-state- 
like wave functions belonging to the continuum spectrum of the Hamiltonian is discussed. Examples 
are shown using realistic nucleon-nucleon forces. 
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I. INTRODUCTION 



Well established methods to treat both, bound and scattering states in A = 3,4 sys- 
tems, are the solution of the Faddeev (A = 3) or Faddeev-Yakubovsky (A = 4) equations 
in configuration or momentum space and the hyperspherical harmonic (HH) expansion in 
conjunction with the Kohn variational principle (KVP). These methods have proven to be 
of great accuracy and they have been tested through different benchmarks [1, 2]. On the 
other hand, other methods are presently used to describe bound states: for example, the 
Green's Function Monte Carlo (GFMC) and No Core Shell Model (NCSM) methods have 
been used in nuclei up to A = 10 and A = 12 respectively [3, 4]. Attempts to use these 
methods in the description of scattering states recently appeared [5, 6]. 

The possibility of employing bound state techniques to describe scattering states has 
always attracted particular attention [7]. Recently, continuum-discretized states obtained 
from the stochastic variational method have been used to study single channel a + n scat- 
tering scattering [8]. The extension to treat coupled channel scattering is given in Ref. [9]. 
In those approaches, the tangent of the phase-shift results in a quotient of two numbers. In 
the former the numerator and denominator are obtained from two integral relations after 
projecting the Schrodinger equation, whereas in the latter the numerator results from an 
integral relation derived by means of the Green's function formalism and the denominator 
from the normalization of the continuum-discretized state. 

Recently two integral relations have been derived from the KVP [10]. It has been shown 
that starting from the KVP, the tangent of the phase-shift can be expressed as a quotient 
where both, the numerator and the denominator, are given as two integral relations. This is 
similar to what was proposed in Ref. [7], however the variational character of the quotient 
and its strict relation with the KVP were not recognized. In fact, it is this characteristic 
that makes possible many different and interesting applications of the integral relations. For 
example, in Ref. [11], the integral relations have been used to compute phase-shifts from 
bound state like functions in the A = 2, 3 systems using semirealistic interactions. Both 
n — d and p — d scattering were considered. The latter process is of particular interest since 
p — d scattering at low energies has been a subject of intense investigations. Initially, the 
Faddeev method has been applied mainly to the neutral n — d reaction. Applications to p — d 
zero-energy scattering were studied in configuration space by the Los Alamos-Iowa group 
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using s-wave potentials [12] and realistic forces [13]. In those calculations the KVP was used 
to correct the first order estimate of the scattering length after solving the Faddeev equations 
in which the partial wave expansion of the Coulomb potential was truncated. Low-energy 
p — d elastic scattering has been studied using the pair correlated hyperspherical harmonic 
(PHH) expansion [14, 15] as well. A benchmark between these two techniques was given 
in Ref. [16]. A different way to treat the Coulomb potential in few-nucleon scattering was 
proposed in Ref. [17], based on the works of Ref. [18], in which the Alt-Grassberger-Sandhas 
equations were solved using a screened Coulomb potential and then the scattering amplitude 
was obtained after a renormalization procedure. A benchmark for elastic p — d scattering 
up to 65 MeV between this technique and the PHH expansion using the KVP has been 
performed [19]. 

Summarizing, the description of scattering states using very accurate methods are at 
present circumscribed to A < 4 systems. On the other hand, accurate methods to describe 
bound states beyond the A = 4 mass system exist. Therefore the discussion of new methods 
to extend these approaches to treat scattering states is of interest. In this discussion the 
treatment of the Coulomb interaction cannot be neglected. In the present work we would 
like to show a detailed application of the integral relations derived from the KVP in which 
A = 3, 4 bound-state-like wave functions are used to compute the scattering matrix using 
realistic nucleon-nucleon NN potentials. In particular, we face the problem of constructing 
degenerate bound state wave functions at a given energy E belonging to the continuum 
spectrum of the Hamiltonian. In fact, in the A = 3 system, the elastic scattering matrix is a 
2x2 matrix for J 71 = l/2 ± and a 3 x 3 matrix for all the other states. This means that, at 
energies below the deuteron breakup threshold, there are two (for J n = l/2 ± ) or three (for 
J > 1/2) scattering states, at the same energy, differing in their asymptotic structure. For 
example, in the J n = l/2 + state, two different asymptotic structures exist corresponding to 
(L, S) = (0, 1/2) or (2, 3/2), being L the relative angular momentum between the deuteron 
and the incoming nucleon and S the total spin 5*. Therefore, at a given energy, the two 
scattering states have a particular combination of the two different asymptotic structures 
determined by the scattering matrix. In the present paper we discuss how to construct 
degenerate bound-state wave functions at a particular energy, belonging to the continuum 
spectrum of the Hamiltonian. Moreover, these states will be used in the integral relations to 
compute the scattering matrix. Examples using realistic forces in the A = 3, 4 systems will 
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be shown. We expect that this study will serve as a guide for calculating scattering states 
in systems with A > 4. 

The paper is organized as follow. In the next section general A = 3,4 bound-state wave 
functions are constructed using the HH expansion. In Section III, a brief derivation of the 
KVP given in terms of the integral relations is discussed. Applications to the A = 3, 4 
systems are shown in Section IV whereas the conclusions are given in the last section. 

II. ,4 = 3,4 BOUND LIKE STATES WITH ARBITRARY J n VALUES 

Following Refs. [14, 20, 22, 23] we give a brief description of a general three- and four- 
nucleon bound state in terms of the hyperspherical harmonic basis. In the case of A = 3, a 
bound-state wave function can be written as a sum of three amplitudes 

* = ^(Xi, Yi) + ^(xj, yj) + ^(xfc, Yk) , (1) 

where Xj,yj are the internal Jacobi coordinates which are defined in terms of the particle 
coordinates as 

^i = Tj-r k , yi = ^( r j + r fe - 2r i) ■ ( 2 ) 

Each z-amplitude has total angular momentum and parity J n and third component of the 
total isospin T z . Using LS coupling, it can be decomposed into channels 

ip(xi,Yi) = ^2(j> a (xi,yi)y a (jk,i) (3) 

a 

y a (jk,i) = {[YeAmLAm A y a %] Sa }^ K%] TaT , (4) 

where Xi,yi are the moduli of the Jacobi coordinates. Each a-channel is labeled by the 
angular momenta £ a , L a , coupled to A a , and by the spin (isospin) {t£) of the pair j, k, 
coupled to the spin (isospin) of the third particle s l a (t l a ) to give S a (T a ). N c is the number of 
channels taken into account in the construction of the wave function and should be increased 
until convergence is reached. The antisymmetrization of the state requires that £ a + s£ + 1 3 ^ 
be odd, while the parity of the state is given by £ a + L a . 

Defining the hyperradius and hyperangle in terms of the moduli of the Jacobi coordinates 

x i = pcos(/) i , y i =psm(/) i (5) 
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the two-dimensional spatial amplitudes can be expanded in terms of the PHH basis as 
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where the hyperspherical polynomials are 
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]\fia,L a - g a norma ii za tion factor, P^" is a Jacobi polynomial and K = l a + L a + 2n is the 
grand orbital quantum number which runs from its minimum value Kq = £ Q + L a to its 
maximum selected value K a . Therefore, the number of hyperradial functions per channel is 
M a = (K a — Ko)/2 + l. The inclusion of the pair correlation function f a {xi) in the expansion 
of Eq. (6) accelerates the convergence taking into account the correlations introduced by the 
strong repulsion of the NN potential (see for example Ref. [20]). 

In the case of the four-nucleon system we use the HH expansion as described in Ref. [21]. 
The wave function having total angular momentum J and parity n can be cast in the form 

* = EE w M^ ] , (8) 

[K] a 

where [K] = K, A, S, T and # are the channel HH-spin-isospin functions having grand 
angular momentum K, orbital angular momentum A, coupled to total spin S, to give a 
total angular momentum JJ Z , and total isospin T. The channel index a labels the possible 
choices of hyperangular, spin and isospin quantum numbers, namely 

a = {£x,£ 2 , 4, L 2 , n 2 , n 3 , S a , S b , T a , T b } , (9) 

compatibles with the given values of K, A, S, T, J and ir. The channel function is 
constructed as a linear combination of the following basis elements 
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Here, 3^V 2 '4 £ 2 re2 na (^) * s ^ ne f our - nuc l eon HH state and Sj (U) denotes the spin (isospin) 
function of particle % and Q indicates the set of the four-nucleon hyperangular variables. 
The total parity of the state is given by ir = (— i)^i+ £ 2+^3_ 

In the present work the A = 3,4 hyperradial functions u? K Ap) are taken as linear combi- 
nations of Laguerre polynomials multiplied by an exponential function: 



u 



[K] 



where A a \K], m are coefficients to be determined and the indeces a, [K] label either a three- 
nucleon or a four-nucleon channel. The polynomials depend on the variable z = ftp, with 
ft a nonlinear variational parameter. Let us define \a, [K],m > as a totally antisymmetric 
element of the expansion basis for the A = 3, 4 systems. In terms of the basis elements, the 
bound-state wave functions given in Eqs.(l) and (8) can be written as 

*„= J2 An a ,[K ]>m \o £ ,[K],m>. (12) 

a,[K],m 

The index n indicates the level of the state with energy E n . The linear coefficients A n a ^ m 
of the wave function and the energy of the state are obtained by solving the following 
generalized eigenvalue problem 

E K, j[K , ]>m ,<a,[K]MH-E n \a',[Klm>>=0 . (13) 

a',[K'],m' 

In the latter equation the dimension of the involved matrices is related to three indices: the 
number of a-channels N c , the number of hyperspherical polynomials for each channel M a 
and, Nl the number of Laguerre polynomials included in the expansion of the hyperradial 
functions of Eq.(ll). The convergence properties of the expansion is analyzed by increasing 
the indices K, m and studying the stability obtained for different values of the nonlinear 
parameter ft. The ground state of the three- nucleon system has total angular momentum 
and parity J w = l/2 + and with N c = 18 an accuracy of 1 keV is reached [22, 24]. The 
corresponding dimension of the PHH basis is D « 2200, considering M a = 8 for the first 
8 channels, M a = 6 for the successive six channels, M a = 4 in the last ones and including 
Nl ~ 20 Laguerre polynomials in the description of the hyperradial functions. After the di- 
agonalization of the whole matrix, D eigenvalues are obtained. The lowest one corresponds 
to the three-nucleon ground state and, with the very extended basis used, it shows a notice- 
able stability with ft. A certain number of negative eigenvalues verifying E n > Ed (with Ed 
the deuteron energy) also appear. Defining the positive energy E® = E n — E d , the corre- 
sponding eigenvectors ty n approximately describe a scattering process at the center-of-mass 
energy E®, though asymptotically they go to zero. The eigenvalues E n present a monotonic 
behavior with ft, as shown in the left panel of Fig. 1, where the AV14 NN potential [25] has 
been used. 

In the right panel of Fig. 1 the lowest eigenvalues obtained from a diagonalization of the 
J n = l/2~ state are shown. As expected this state is not bound, though several negative 

6 



states appear with energies in the interval Ed < E n < 0, characterized with a monotonic 
behavior with 0. As before, these states approximately describe a scattering process at the 
center-of-mass energy E®. In Fig. 1 the deuteron energy is indicated by the dotted-dashed 
line whereas the three dashed lines indicate the lab energies Ei^ = 1, 2, 3 MeV, respectively. 
Interestingly, the energies of the J n = 1/2" state appear in pairs. This can be understood 
noticing that the J 71 = 1/2" scattering states are twofold degenerate at a given energy, as 
the scattering matrix has dimension of two. This degeneration arises from the two possible 
asymptotic configurations in which the relative angular momentum of the deuteron and the 
third nucleon is L — 1 and the total spin can take the values 5 = 1/2 and 3/2. Also the 
J 71 " = l/2 + state is twofold degenerate, having two possible asymptotic configuration with 
the values L — 0,S — 1/2 and L = 2,S = 3/2. However, in this case, the different L 
values produce different contributions to the kinetic energy with the consequence that the 
two degenerate states appears with a larger separation compared to the J n = 1/2" case. 
However, this difference reduces as the basis is enlarged. 

To analyze further the hypotesis that the states organize in pairs corresponding to the 
two different asymptotic configurations in both J 71 = 1/2^ states, in Table I the different 
occupation probabilities are given. For the J w = l/2 + state the occupation probabilities of 
the S- P and P-waves, P$, Pp and Pd, have been computed. The Eq level corresponds to 
the ground state and the successive levels organize in mostly S-wave (Ei and E 3 ) and mostly 
P-wave (E2 and E4) states, alternatively. In the case of the J 77 = 1/2" state, the occupation 
probabilities of the P-wave with total spin values S = 1/2 and 3/2, Pp 2 and Pp 7 ' 2 , as well 
as Pd have been computed. From the table we can observe that the levels organize in pairs, 
being one of the states mostly a P-wave state with 5" = 1/2 and the other mostly a P-wave 
state with S = 3/2. This organization is indicated in Fig. 1 with colors. For the J n = l/2 + 
the E level, shown as a black solid line, is practically constant with 0. The levels with high 
L = (L = 2) occupation probability are given in red (blue) respectively. In the case of the 

3 /2 1 /2 

J w = 1/2", the levels with high P p (P p ) probabilities are given in red (blue) respectively. 
For small values of /3 the spectrum tends to be denser since, in this case, the polynomials 
can contain more oscillations before the action of the exponential tail becomes significant. 
As /3 increases the number of negative eigenvalues decreases. In the case in which a bound 
state exists, as in the case of the J 71 " = l/2 + state, the basis is sufficiently large to guarantee 
a correct description of it as the control parameter /3 is varied. As we will see, the wave 
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functions \l/ n corresponding to energy levels Ed < E n < 0, can be used to determine the 
scattering matrix at specific energies. 

In the case of the A = 4 system we analyze the single channel J n = + state with 
T — T z — 1, corresponding to the p — 3 He system. Using the N3LO-Idaho potential [26], 
the Hamiltonian matrix has a total dimension D »s 84000, obtained expanding the wave 
function on the HH basis, as as previously described, with K max = 44, corresponding to 
about 3500 HH states, and Nl = 24. For this values of D, the matrix can be diagonalized 
using standard iterative methods. In Fig. 2 the first eigenvalue is shown as a function of 
the control parameter /3. Clearly the lowest eigenvalue is above the 3 He threshold, fixed for 
the N3LO-Idaho potential at -7.128 MeV, since four nucleons in the isospin channel T = 1 
does not present a bound state. The three dashed lines correspond to three lab energies 
(3.13, 4.05 and 5.54 MeV) at which experimental data exist. Similar to the previous cases 
in A = 3, we will use these four-body bound state wave functions to determine the p — 3 He 
scattering matrix at the indicated energies. 

III. THE KVP IN TERMS OF INTEGRAL RELATIONS 

Following Refs. [15, 23] a general scattering state with A = 3, 4 can be written as a sum 
of two terms 



The first term, ^c, describes the system when the A nucleons are close to each other. 
For large interparticle separations and energies below the breakup threshold in more than 
two pieces it goes to zero, whereas for higher energies it must reproduce a three or four 
outgoing particle state. It can be written as a sum of amplitudes corresponding to the 
cyclic permutations of the Jacobi coordinates. Each amplitude ^^({xj}) has total angular 
momentum and parity J n and third component of the total isospin T z (here {x^} represents 
the set of Jacobi coordinates with ordering of the particles % for the A = 3 or A = 4 systems). 
For energies below the breakup threshold in three pieces, it can be expanded in terms of the 
totally antisymmetric states 



(14) 




(15) 



a,[K],m 
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The second term, ^a, in the scattering wave function of Eq.(14) describes the relative 
motion of the two clusters in the asymptotic region. For A = 3, ^>a describes the relative 
motion between the deuteron and the incident nucleon, whereas for A = 4 we will limited 
the description to an incident nucleon on 3 He or 3 H. It can be written as a sum of amplitudes 
whose generic form for A = 3 is given by 

fi£ 5J (x J? y 4 ) = ™MK X L { V i) { [[Y la (x t )si k U] s Y L (y t )} Jj2 [^¥]tt, , (16) 

l a =0,2 

where Wi a (xi) is the l a — 0, 2 deuteron wave function, s£ = 1, = 0, and L is the relative 
angular momentum of the deuteron and the incident nucleon. The superscript A indicates 
the regular (A = R) or the irregular (A = I) solution of the Schrodinger equation in the 
asymptotic region. In the p — d (n — d) case, the functions 7Z X are related to the regular or 
irregular Coulomb (spherical Bessel) functions. The functions Q x can be combined to form 
a general asymptotic state 



^L5j( x «' y») + 5J j r ^LL'^L , s"j( x «' y« 



i=l,3 L L'S 



(17) 



where 

fi L5j( x i>yi) = «oo^L5j( x ^yi) +uoitt I LS j(xi,yi) , (18) 
fi L5j( x i ; yi)= wio^f S j(xi,yi) + unfi£ SJ (x i ,y i ) . (19) 

The matrix elements w^- form a matrix tt that can be selected according to the four different 
choices of the matrix C = /C-matrix, /C~ 1 -matrix, S'-matrix or T-matrix. It should be noticed 
that the irregular solution has been opportunely regularized at the origin 

K I L (y) = (l-e™') L+1 G L (y) (20) 

where r^d = (y/3/2) y is the nucleon-deuteron separation and the parameter 7 is fixed 
requiring that TV L {y) = Gl(v) asymptotically. Moreover, Gl(v) is the irregular Bessel 
function or the irregular Coulomb function in the case of n—d or p— d scattering, respectively. 
The description for A = 4 can be found in Ref. [27] 

A general three- or four-nucleon scattering wave function for an incident state with rela- 
tive orbital angular momentum L, spin S, total angular momentum J and energy below the 
three-particle breakup threshold is 

\*1sj>= E A %x\J a >i K \> m> +\ n isj>> ( 21 ) 

a,[K],m 
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and its complex conjugate is ^2sj- ^ variational estimate of the trial parameters in the 
wave function ^^ SJ can be obtained by requiring, in accordance with the generalized KVP, 
that the functional 

[ J C S £] = J Cf v - ^{^IsjW - ElV+yj) , (22) 

be stationary. Applications of the complex KVP for N — d scattering can be found for 
example in Refs. [15, 22, 28]. In the case in which the variational principle is formulated in 
terms of the /C- matrix, we get: 

[ J /Cff,] = J /Cff; - i^l SJ \H - E\V+g,j) . (23) 

Calling the set of indeces //{= a, [K],m} and % = {L, S, J}, the variation of the functional 
[ J K-ll\ — [ fca] w ith respect to the linear parameters A 1 ^ leads to the following two sets of 
linear equations 

< n\H - E\fi' > A°J = -<fi\H- > (24) 
^2<H\H- E\n' > Atf = - < n\H - E\Q\ > , (25) 

in accordance of the two possible asymptotic scattering states and Q}. From the above 
equations the two sets of coefficients A ' 1 , A 1 ' 1 can be obtained. Furthermore, multiplying 
the sets by these coefficients and summing on /x, it is possible to reconstruct the scattering 
state and the above equations can be formally cast as 

< V C \H-E\#t >=0. (26) 

The variation of the functional with respect to the linear parameters /Qj results 

Sij- < Q]\H - E\y+ > - < ®r\H - E\Vt) >= . (27) 

Using the normalization condition 

< n° t \H - E\n) > - < n)\H - e\s% >= 5 tj , (28) 

the scattering wave function verifies 

< V7\H - E\Q) > - < Vt)\H - E\V+ >= 8^ (29) 
< S$\H - E\^f+ >-<Vj\H- E\n° >= Kij , (30) 
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allowing to reduce Eq.(27) to 

< Stij\H-E\V+ >=0. (31) 

The second order estimates of the /C-matrix elements [ J /Cff/] = [ /Cj,-/] are obtained 
replacing in the functional of Eq.(23), the first order solutions given by Eqs.(26) and (31). 
It results 

[ K iV ] = K iV - <^\H- E\V+ > (32) 
that can be further reduced using Eq.(30) to 

[ £..,] = - < yr\H - E\S%, > . (33) 

This final form of the KVP is a direct consequence of the particular form selected for the 
asymptotic scattering state given in Eq.(17) in which the flux of the regular wave Vt® has been 
set to one. As we will see in the following, it is useful to define an asymptotic scattering 
state with general coefficients in both the regular and irregular waves. Accordingly, the 
asymptotic scattering state now reads 



8=1,3 



J2 JA LL> n L>s>A^yi) + Y< JB LL>n 1 LISIJ (x l ,y l ) • (34) 

MS' US' 

The coefficients ^4ff/ and J Bf i f i , form the matrices A and B respectively and the scattering 
matrix results K. = A~ X B. Starting with a scattering state that has this asymptotic behavior, 
the relations of Eqs.(29) and (30) result 

< *r\H - E\n) >-<n]\H- E\$t >= Aj 

< n^H - E\^+ >-<Vj\H- Eltf >= Bij , (35) 

and, using Eqs.(31)-(33), the KVP takes the particular form 

Aij = < Vr\H-E\n) > 
[By] = -<*T\H-E\rtj> 
[JC] = A- l \B] (36) 

where [B] and [JC] are second order estimates. 

Eqs.(36) formulate the KVP in terms of integral relations depending on the internal 
structure of the scattering wave function In fact (H — E)Qj and (H — E)£Vj go to 
zero as each of the three Jacobi coordinates jjk goes to oo (k = 1,2,3), since Q®' 1 are the 
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solutions of (H — E) in that limit. Therefore, in Eqs.(36), it would be possible to use trial 
wave functions that are solutions of (if — E) in the interaction region but do not have 
the physical asymptotic behavior indicated in Eq.(34). In particular, it would be possible 
to use the bound-state wave functions \l/ n described in the previous section to calculate the 
scattering matrix corresponding to a center-of-mass energy E®. This is discussed in the next 
section. 

IV. SCATTERING MATRIX FROM BOUND-STATE-LIKE WAVE FUNCTIONS 

In section II the construction of A = 3, 4 bound states having general quantum numbers 
jit corresponding to different levels with negative eigenvalues E n has been discussed. In 
the case of the A = 3 state J = l/2 + the E level and the corresponding wave function 
describe the energy and structure of the triton or 3 He for the two possible values of 
T z = 1/2 or —1/2, respectively. Using the nonlinear parameter (3 as a control parameter it 
was possible to construct states with eigenvalues E n in the region E^ < E n < 0. In a similar 
way, it is possible to contruct these kind of states for arbiratry values of J n . As an example, 
in section II, the case J = l/2~ was explicitly discussed. Furthermore, it was shown that 
these states organize sequentially having occupation probabilities that can be connected 
with the different components of a scattering state, corresponding to the different values of 
the quantum numbers L, S, J. The number of these components fixes the dimension of the 
scattering matrix and, correspondingly, the degeneration of the state. Therefore, in order to 
construct a scattering state using bound-state-like functions, those components can be taken 
into account considering sequential solutions having the same energy. To this end the control 
parameter /3 can be used to select sequential solutions at the same eigenvalue E n . This is 
shown in Fig. 3 for three different cases. The three dashed lines in both panels of the figure 
indicate the energies corresponding to incident energies in the lab system Ei a b = 1,2,3 
MeV. As explained in section II, the red and blue lines show the variation of sequential 
eigenvalues as a function of (3 having the different structures given in Table I. The circles 
in Fig. 3 indicate the points in wich the eigenvalues cross the dashed lines and, accordingly, 
at those specific values of two different solutions, ^\ and can be constructed having 
the same energy E n and presenting a very different internal structure. These two states are 
solutions of (H — E n )fy^ 2 = in the internal region and, since they are square integrable 
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states, they go to zero asymptotically. However the integral relations of Eq.(36) depend on 
the internal part of the wave function and therefore it would be possible to use and 
as trial wave functions. In this case the second order estimate of the scattering matrix is 

Aij = <%\H-E n \n} > 
[Bij] = -<V n \H-E n \tt° 3 > 
[K] = A- l [B] (37) 

where i,j indicate either the two solutions, ^l; 2 , and the two possible values of the set of 
quantum numbers (L, S, J) in J = l/2 ± . 

For the A = 4 case we have analyzed the single channel J n = + state with T = T z = 1. 
In Fig. 4 we show the three cases (indicated with circles) at which, for specific values of the 
control parameter /3, the eigenvalue matches the selected energies. Accordingly the second 
order estimate of the scattering matrix can be obtained in each 

A = < V n \H- E n \tt l > 
[B] = - < V n \H - E n \n° > 

[K] = A^IB] (38) 

In this case we are considering a single channel state and therefore the scattering matrix 
results a scalar. 

In the following, results of phase-shifts and mixing parameters for the n — d system, 
calculated using the AV14 NN potential, are presented for the J w = l/2 + state in Fig. 5, 
and for the J n = l/2~ state in Fig. 6, at the three selected energies E^ = 1, 2, 3 MeV. The 
stability of the results with 7, the regularization parameter introduced in Eq.(20), is chosen 
as a convergence criterion. This criterion has been discussed in Refs. [10, 11] and essentially 
it establishes the quality of as solution of (H — E n )W n = 0. In fact, if is a good 
solution, the integrals of Eq.(37) are largely independent of 7. The results are compared to 
the benchmark of Ref. [1] given in the figures as a red line. The results of the application 
of Eq.(37) are shown as filled circles corresponding to values of 7 varying from 0.25 fm -1 to 
1.25 fm _1 . We can observed a good stability on this interval and, furthermore, the results 
are in very good agreement with those of Ref [1]. 

In Fig. 7 results are given for the n — dJ^ = 3/2 + state. In this case the /C- matrix is a 3 x 3 
matrix, corresponding to asymptotic configurations having L = 0, S = 3/2, L — 2, S — 1/2 
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and L = 2, S = 3/2. The diagonalization of the Hamiltonian matrix in the J 71 = 3/2 + case 
produces sequential eigenvalues with occupation probabilities in accordance with these three 
configuration. Using the control parameter f3, three sequential eigenvalues can be chosen to 
have a particular E n value as has been done before for the J n = 1/2^ states. In the figure 
we observe a good stability with the regularization parameter 7 and a close agreement with 
the results of the Ref. [1]. In Fig. 8 results for the p — d J = l/2 + state are given. The 
description of the p — d process at low energies presents some problems using the Faddeev 
equations. In Ref. [16] a benchmark for p — d scattering has been produced using the HH 
method and the Faddeev method in configuration space. The results of the benchmark are 
shown as a red line in Fig. 8. From the figure we can observe that the results using the 
integral relations reproduce extremely well the benchmark results. This is an important 
point since in bound-state-type calculations the treatment of the Coulomb potential does 
not present any troubles. 

The p — 3 He results are given in Fig. 9 for the three selected energies. The phase-shift for 
the + state is shown as a function of the regularization parameter 7 (filled circles). As a 
comparison, the results of the recent benchmark of Ref. [29] are shown as a red line. We can 
observe a good stability with 7 indicating that the four-nucleon bound-state eigenfunction \1> 
is a good solution of (H — E)^ = at the specified energies. Moreover the results obtained 
using the integral relations are in close agreement with those of the benchmark. 

V. CONCLUSIONS 

In this work the elastic scattering matrix has been determined using bound-state-like 
wave functions. To this end two integral relations derived from the KVP have been used. 
Initially, these integral relations were derived in Ref. [10] in order to extract phase shifts 
from the solutions calculated using the hyperspherical adiabatic expansion in the three- 
nucleon system. In this method the boundary conditions at large distances are imposed in 
terms of the hyperradius p. However, as explained in Ref. [10], the boundary conditions 
depend explicitly on the Jacobi coordinates Xj, y, describing the asymptotic configuration of 
a deuteron formed by particles (J, k) and an incoming nucleon (particle i). The equivalence 
between imposing the boundary conditions in p or in the Jacobi coordinates directly, results 
at very large values of the hyperradius where the relation p m is verified. As a consequence, 
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the phase shifts obtained from the adiabatic expansion require a large number of terms to 
converge. On the other hand, the phase shifts obtained as a quotient of the two integral 
relations converge much faster and, in fact, the rate of convergence is similar to that one 
obtained in the case of bound state solutions. The reason behind this fact is that the 
integral relations depend only on the internal part of the wave function. Therefore, it is 
enough that the wave function verifies (H — E)^ = in the internal region to obtain almost 
exact results for the scattering matrix at the center-of-mass energy E. This characteristic 
allows to apply the integral relations using bound-state-like wave functions obtained from 
a direct diagonalization of the Hamiltonian H . Eigenvectors corresponding to eigenvalues 
belonging to the continuum spectrum of H can be used as inputs to determine the scattering 
matrix at fixed values of E. Applications for single-channel solutions using semi-realistic 
NN potentials are given in Ref. [11]. The coupled- channel case of an atom colliding a dimer 
formed by other two atoms is given in Ref. [30]. 

In the present work, applications to elastic scattering of a nucleon on a deuteron {A = 3) 
or on 3 He (A = 4) below the breakup threshold, using realistic nucleon-nucleon potentials 
has been discussed. In particular, for A = 3, two or three solutions at the same energy 
have to be determined corresponding to the different possible asymptotic configurations of 
the system. A detailed construction of such solutions, using the nonlinear parameter /3 as 
a control parameter, has been analyzed. Moreover it has been shown that the eigenvectors 
of successive eigenvalues organize in pairs (for J = l/2 ± ) or in triplets (for J > 1/2), 
corresponding to the different asymptotic structures. The control parameter (3 has been 
tuned to find solutions having the same eigenvalue that have been usend to calculate the 
scattering matrix at the selected energy. The obtained results are in close agreement with 
those presented in the A = 3 benchmarks of Refs. [1, 16] and in the A = 4 benchmark of 
Ref. [29]. In particular, the results for p — d and p — 3 He scattering demonstrate that the 
scattering matrix can be calculated using bound-state-like wave functions also in scattering 
of charged particles. 

Well established bound-state methods to diagonalize the nuclear Hamiltonian in systems 
with A > 4 already exist. The formulation of the scattering matrix presented in this work 
will allow to extent these studies to the low-energy continuum spectrum. It will be then 
possible to compare theoretical predictions for scattering observables to the experimental 
data, in order to evaluate the capability of the present models for the interaction to describe 
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the nuclear structure. Studies along this line are at present intensively pursued. 
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J 


= 1/2+ 


Ps{%) 


Pp(%) 


Pd{%) 




Eq 


90.96 


0.08 


8.97 




Ei 


90.04 


0.00 


5.96 




E2 


1.22 


2.72 


96.06 




E3 


94.20 


0.00 


5.80 




E A 


1.21 


2.70 


96.09 


J 


= 1/2" 


Pp /2 (%) 


P 3 p /2 (%) Pd{%) 




Ei 


3.36 


94.87 


1.77 




E2 


93.80 


3.42 


2.78 








94.71 


1.82 




E A 


93.42 


3.86 


2.72 




E 5 


4.32 


93.87 


2.01 




Ec, 


92.72 


4.77 


2.51 



TABLE I: Occupation probabilities of the different states shown in Fig. 1 at f3 = 2.5 fm 1 (J 71 " = 
1/2+) and (3 = 2.0 fm" 1 (J 71 = 1/2 - ). 
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FIG. 1: (Color online) The lowest A = 3 eigenvalues, using the AV14 potential, for the J w = l/2 + 
(left panel) and J n = l/2~ (right panel) states, as a function of the nonlinear parameter /3. 
The solid black line (left panel) represents the triton energy whereas the colored lines indicate the 
eigenvalues embedded in the continuum as explained in the text. The dotted-dashed line represents 
the deuteron energy, the dashed lines indicate the labs energies Ei a b = 1,2,3 MeV. 
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FIG. 2: (Color online) The lowest A = 4 eigenvalue, using the N3LO potential, for the J n = + 
state with T = T z = 1 as a function of the nonlinear parameter /3. The colored lines indicate the 
eigenvalues embedded in the continuum as explained in the text. The dotted-dashed line represents 
the triton energy, the dashed lines indicate the labs energies E[ a ^ = 3.13,4.05,5.54 MeV. 
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FIG. 3: (Color online) The same as Fig.l in which two sequential solutions having the same 
eigenvalue are selected (indicated by the circles) in the three cases corresponding to incident energies 
Eiab = 1, 2, 3 MeV respectively. 
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FIG. 4: (Color online) The same as Fig. 2 in which the lowest eigenvalue is selected (indicated 
by the circles) in the three cases corresponding to the lab energies E cm = 3.13,4.05,5.54 MeV 
respectively. 
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FIG. 5: (Color online) The n — d J w = l/2 + phase-shifts and mixing parameters as a function of 
the regularization parameter 7 at the three indicated energies. The red line corresponds to the 
results of Ref. [1]. 
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FIG. 6: (Color online) The n — d J = l/2~ phase-shifts and mixing parameters, for the AV14 
potential, as a function of the regularization parameter 7 at the three indicated energies. The red 
line corresponds to the results of Ref. [1]. 



24 



R =1 MeV 

lab 



R = 2 MeV 

lab 



R =3 MeV 

lab 



r-, 0.5 
T3 



i° -0.5 

r, 



-1 

-1.5 



a -47 

+ 



-a 



■47.5 
1 



a> 0.5 

C3 





-0.5 



1 1 1 1 1 














— 4 

1 D 3/2 






i 1 i 1 i 


I 

2 s 

- a l/2 


i 


i 


i 1 i 1 i 


i 

_ ^3/2+ 


i 


i 








— • — 


; ; 


_ £ 3/2+ 






^3/2+ | | 
i 1 i 1 i 



0.5 
y [fm" 1 ] 



1 



-2 
-4 

61.2 

■61.4 

61.6 
1 

0.5 


-0.5 



- ' 2 ' ' 

- D 3/2 


1 


4 

- D 3/2 




1 1 


1 


1 1 1 1 

" 4 s 

^3/2 


1 






1 1 

^3/2+ 


1 




~~ £ 3/2+ 




" ^3/2+ 




■ 1,1, 



1.5 



0.5 



1 



2 



-2 

-4 
-70 

■70.5 

-71 
1 

0.5 


-0.5 



"D 



3/2 



D 



3/2 









— 1 — 
" 4 s 

^3/2 


1 


1 








_ s 3/2+ 







•'3/2+ 



11 



3/2+ 

T 1 T 1- 



1.5 



0.5 



1 



1.5 



Y [fm" ] 



y [fm" ] 



FIG. 7: (Color online) The n — d J = 3/2 + phase-shifts and mixing parameters, for the AV14 
potential, as a function of the regularization parameter 7 at the three indicated energies. The red 
line corresponds to the results of Ref. [1]. 
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FIG. 8: (Color online) The p — d J = l/2 + phase-shifts and mixing parameters, for the AV14 
potential, as a function of the regularization parameter 7 at the three indicated energies. The red 
line corresponds to the results of Ref. [16]. 
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FIG. 9: (Color online) Thep- 3 He J = 0+ phase-shift, for the N3LO potential, as a function of the 
regularization parameter 7 at the three indicated energies. The red line corresponds to the results 
of Ref. [29]. 
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